Electronic Structure of Sr 2 FeMoO 
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We have analysed the unusual electronic structure of Sr2FeMoC>6 combining ab-initio and model 
Hamiltonian approaches. Our results indicate that there are strong enhancements of the intraatomic 
exchange strength at the Mo site as well as the antiferromagnetic coupling strength between Fe and 
Mo sites. We discuss the possibility of a negative effective Coulomb correlation strength {U e ff) at 
the Mo site due to these renormalised interaction strengths. 



Colossal magnetoresistance in transition metal oxides 
[1| has attracted a great deal of attention in recent times 
[2| owing to their technological as well as fundamental 
importance. The systems based on doped manganites 
have been most exhaustively studied and their behaviors 
have been broadly understood in terms of the double- 
exchange mechanism |^ , possibly modified by the pres- 
ence of various kinds of polarons j5j. Recently, an or- 
dered double perovskite oxide without any manganese, 
Sr2FeMo06, has been shown Q to have CMR at a low- 
field even at room temperature. This spectacular MR in 
Sr2FeMoC>6 is intimately connected with a large magnetic 
T c (?« 450 K) of this compound. The essential physics of 
Sr^FeMoOg is believed to be very similar to that of the 
manganites. However this apparent similarity between 
the manganites and S^FeMoOg is surprising, if not to- 
tally unexpected. 

Mo is not a strongly correlated system and conse- 
quently, a magnetic moment at the Mo site is a rar- 
ity, in contrast to Mn. In the ordered double perovskite 
Sr2FeMo06, the transition metal sublattice sites are oc- 
cupied alternately by Fe and Mo ions. Experimental 
evidence of ferrimagnetic order as well as band struc- 
ture results establish an antiferromagnetic coupling be- 
tween Mo and Fe atoms, in contrast to the ferromagnetic 
coupling in manganites. The large magnetic transition 
temperature in S^FeMoOg, however, points to a large 
interatomic exchange coupling strength, J, between Fe 
and Mo, comparable to that between the Mn-Mn pairs 
in the manganites, inspite of the normally expected non- 
magnetic nature of Mo. Mo is usually nonmagnetic as 
the intraatomic exchange strength, /, within the Mo 4d 
manifold is small, typically 0.1-0.2 eV, nearly an order of 
magnitude smaller than that in the 3d transition metal 
systems |?]] . Additionally, the Mo Ad bandwidth tends to 
be substantially larger compared to the exchange split- 
ting. However, the magnetic structure of this compound 
requires a strong spin-splitting of the delocalized bands, 
leading to nearly complete spin polarization. This is- 
sue is moot in the present context, since the spectacular 
properties of this system - a high temperature and low- 
field MR - are critically dependent on the large T c and 
the half-metallic ferromagnetic state. We investigate the 



origin of these apparent anomalies by combining state-of- 
the-art ab initio methods as well as realistic multiband 
many-body models for electronic structure calculations. 
Our results clearly establish an unusual renormalization 
of the intraatomic exchange strength at the Mo sites, 
arising from Fe-Mo interactions; this causes the effective 
I to be much larger than the expected atomic value. We 
further evaluate J to show that it is indeed large, ow- 
ing to the renormalization of /. Finally, we suggest the 
novel possibility of a negative effective Coulomb inter- 
action strength, U e ff (= U - I), for the Mo sites due 
to the enhancement of /, with some preliminary experi- 
mental evidence for the same. Our results underline the 
unusual and unique aspects of the electronic structure of 
Sr 2 FeMo06, that have not been appreciated so far. 

We have computed the band dispersion of geometry 
optimised |6| S^FeMoOg within the framework of linear 
muffin-tin orbital (LMTO) method || using the gener- 
alized gradient approximation (GGA) for the exchange- 
correlation part. Figs. 1(a) and 1(b) show the orbital 
projected band dispersions of S^FeMoOg for the major- 
ity and minority spin channels; here the fatness of the 
bands in each panel is the weight of the indicated orbital 
in the wavefunction. The minority down spin bands cross 
the Fermi level, while the majority up spin bands exhibit 
a bandgap of ~ 0.8 eV. The bands below -2 eV are pre- 
dominantly of oxygen character while bands crossing the 
Fermi level and ranging from -2 eV to about 2 eV have 
significant mixing between the Fe-ci and Mo-c? characters 
with some small admixture of oxygen p states. The pres- 
ence of approximate cubic symmetry of the octahedral 
co-ordination of the oxygen atoms around the transition 
metal sites, results in a splitting of the d levels into t^g 
and e g orbitals. For the minority spin channel, the Fe 
t2 g and Mo t2 g bands are partially filled while Fe e g and 
Mo e g bands remain empty. The occupied part of the 
bands near the Fermi level in majority spin channel [see 
Fig. 1(a)] are mainly composed of Fe d character which 
hybridizes with the oxygen p states. The narrow bands 
lying immediately above the Fermi level spanning an en- 
ergy range of about 0.7 eV to 1.2 eV are predominantly 
from Mo t2 g contribution, while the Mo e g bands are 
further high up in energy. In order to estimate the elec- 
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tronic interactions strengths, we fitted the ab initio band 
dispersions, shown in Fig. 1, in terms of a tight-binding 
model containing d orbitals at the Fe and the Mo sites 
and the p orbitals at the oxygen sites. Such methods have 
been successful in obtaining realistic estimates of the in- 
teraction strengths including that for the intraatomic ex- 
change interactions 0,^] . This analysis yields a value of 
the intraatomic exchange splitting strength at the Mo site 
to be 0.13 eV, estimated from the spin polarization of the 
Mo site energies; this estimate is consistent with our ex- 
pectation of a small / at the Mo sites, but is in apparent 
contradiction to the enhanced T c and strong polarization 
of the bands in S^FeMoOe. It is intriguing at this stage 
to note that the up-spin bands related to the Mo d states 
appear in the energy window centred at ~ 1 eV, while 
the down-spin states appear near the Fermi level with a 
larger bandwidth suggesting a spin-splitting of the Mo d 
band, considerably larger than the bare /. This clearly 
indicates that the bare I at the Mo site must be strongly 
renormalised giving rise to an enhanced I e // leading to 
the pronounced spin polarization of the Mo d bands. In 
order to understand this renormalisation process, as well 
as to estimate the effective I in contrast to the bare / 
at the Mo sites, it is necessary to integrate out the other 
degrees of freedom, such as the Fe d states, from the ab 
initio results, so that the few orbital description retains 
the effect of renormalisations as a change in the effective 
interaction strengths. 

Very recently, the third generation LMTO method |l(]] 
has been proved to be a powerful tool for deriving few 
orbital, orthogonal tight-binding Hamiltonians starting 
from the all orbital Hamiltonian representation used for 
the local density approximation or GGA self-consistency 
procedure, by integrating out other degrees of freedom. 

In order to estimate the I e // at the Mo sites, responsi- 
ble for the large spin polarization of the bands, we have 
used this technique to derive the effective Mo tig Hamil- 
tonian from the all orbital Hamiltonian. Figs. 2(a) and 
2(b) show the energy bands in up and down spin channels 
respectively, obtained by diagonalising the 3x3 Hamilto- 
nian defined in the effective Mo ti g basis. It is clear that 
this very powerful approach captures the dispersions of 
the Mo tig orbitals very effectively. The spin-splitting of 
the on-site Mo tig energies can be directly read off from 
the diagonal terms of the real space representation of the 
Hamiltonian matrix (obtained by Fourier transformation 
of the downfolded Hamiltonian in k-space H(k) — > H(R)) 
and was found to be 0.8 eV. For fully polarized Mo ti g 
electrons this is then the ab-initio estimate of I e ff, while 
taking into account the LMTO estimate for the magne- 
tization (Nj-Nj,=0.34 /is) increases the 7 e // value to ~ 
2.3 eV. This spectacular enhancement of the I e ff com- 
pared to the bare I of Mo is in fact easy to rationalize in 
terms of the details of the electronic interactions in such 
a double perovskite system. Within the band structure 
calculations performed for the nonmagnetic ground state 



we find the Mo 4<i states to be located at higher energies 
compared to the Fe 3c? states. By fitting the nonmag- 
netic band structure to a tight-binding model, this en- 
ergy separation is found to be about 1.4 eV. Considering 
the ferrimagnetic arrangement of the Fe and Mo sites, the 
orbital energies are as shown in the left panel of Fig. 3 
in the absence of any hopping interaction. The Fe 3+ site 
has a larger exchange splitting compared to the crystal 
field splitting, while the situation is reversed for the Mo 
site. In the presence of hopping interactions, there is fi- 
nite coupling between states of the same symmetry at the 
Fe and the Mo sites, leading to perturbations of the bare 
energy levels. It is then easily seen that the Mo ti g ^ state 
will be pushed up and the Mo ti g state will be pushed 
further down by hybridization with the corresponding Fe 
states, as shown in the figure. These opposite movements 
of the Mo up- and down-states increase the energy sep- 
aration between these two states, thereby substantially 
increasing the effective exchange splitting at the Mo site. 
Thus, the renormalisation of I e ff at the Mo site is driven 
by the large I at the Fe site and the substantial hopping 
interaction coupling the two sites. It is interesting to 
note, and also follows from the arguments given above, 
that the antiferromagnetic coupling between Fe and Mo 
is crucial for the enhanced I e tf and consequently, the T c . 
A ferromagnetic arrangement on the other hand would 
result in a reduced T c , 

The above analysis as well as the estimate of I e // are 
based on effectively single-particle theories. In order to 
verify that the present conclusions are not seriously af- 
fected by well-known limitations of such theories, we have 
also used many-body calculations to estimate I e ff- In 
order to obtain a realistic description of the specific sys- 
tem, we fixed the hopping interactions between the var- 
ious orbitals to those estimated from the tight-binding 
fit of the full ab initio calculations. Furthermore, we 
performed Hartree-Fock calculations for the multiband 
model for the lattice involving Fe and Mo d and oxy- 
gen p states and reproduced the half-metallic ferromag- 
netic state in order to obtain appropriate estimates for 
the charge transfer energies, intraatomic Coulomb and 
bare exchange strengths. The parameters used are 0.65, 
1.34, -2.33, and 1.47 eV for oxygen ppa, oxygen-Mo pdir, 
oxygen- Fe pda and pdn hopping strengths; £/p e , Ip e , 
Umo, Imo, Ajr e , and Aj,(„ were estimated to be 5.0, 0.8, 
1.0, 0.2, 2.8 and 4.0 eV, respectively. Additional crystal 
field splittings between the ti g and e g orbitals at the Fe 
site and Mo site were introduced. 

For calculation of I e ff, we consider a multiband 
Hubbard-like Hamiltonian for the finite cluster Fe-O-Mo- 
O-Fe. The Hamiltonian includes on-site energies, hop- 
ping interactions and the on-site Coulomb interactions 
at Fe and Mo sites. The full Hilbert space for even such 
a small cluster has a dimension of ~ 2xl0 9 ; since such 
a large calculation is not tractable, we included the hop- 
ping interactions connecting only the transition metal 



2 



t 2g orbitals to the oxygen p^ orbitals. This is a rea- 
sonable approximation, as the relevant charge dynamics 
involves only these orbitals with the e g orbitals partic- 
ipating in Coulomb interactions. Moreover, we restrict 
ourselves to the S z =9/2 subspace in conformity with 
the experimentally observed antiferromagnetic arrange- 
ment of the transition metal moments. The renormalizcd 
I e ff can be estimated from this cluster by calculating 
5E/8n^ — SE/Sn^, where E is the total energy and «j 
and ni are the corresponding Mo occupancies. For the 
above set of parameter values, the I e f / is calculated to 
be 0.95 eV with integral charge (i.e. <$n=l) fluctuations, 
exhibiting a large enhancement from the bare / of 0.2 
eV used in the calculation and in very good agreement 
with the estimate obtained from the ab initio results. 
Since this estimate depends on the choice of the param- 
eter strengths, we have calculated I e ff with reasonable 
variations of the input parameters; we find that I e ff is 
always between 0.8 and 1.5 eV. 

In order to calculate the J, coupling the Fe and the Mo 
sites antiferromagnetically, we perform Hartree-Fock cal- 
culations with the transition metal d-oxygen p multiband 
model with the parameter values specified earlier. We 
considered different magnetic ground states correspond- 
ing to different values of the spin-density wave vector q. 
Small rotations of the spin moments about the ground 
state gave us J(q) plj. We then mapped the results on 
to a nearest neighbor Heisenberg model and obtained a 
value of 18 meV for the coupling between Fe and Mo. 
A similar analysis [[l2| for the half-metallic 30% hole- 
doped LaMnC>3 gave us 12.5 meV. Most importantly, we 
find that the J in S^FeMoOg is in fact somewhat larger 
than that in the manganites, explaining the enhanced 
magnetic transition temperature in the ordered double 
perovskite system compared to the hole doped mangan- 
ites (< 300 K). This enhanced value of J is of course 
intimately connected with the strong spin-polarisation of 
the Mo d states via the renormalised I e ff at the Mo sites. 

It is to be noted that the exchange strength I e f f at the 
Mo site is comparable or larger than the expected value 
of the direct Coulomb interaction strength, U, at the Mo 
site. This leads to the intriguing situation that the ef- 
fective electron-electron repulsion strength, U e ff = U - 
I e ff, for Mo may become negative. If this were indeed to 
be true, there should be some measurable consequences 
in terms of the electronic properties. In order to verify 
this, we have performed core level x-ray photoemission 
study of the Mo 3d core level in S^FeMoOg. The core 
level spectrum is shown in Fig. 4. For comparison, we 
also show the Mo 3d spectra from two reference mate- 
rials, Mo metal and M0O3 in the insets A and B. The 
insets clearly show that the core level spectra from refer- 
ence Mo systems contain only a doublet feature arising 
from the spin-orbit splitting of the core level into 3d 5 / 2 
and 3<i 3 /2 states. In sharp contrast, the Mo 3d spectrum 
from Sr2FeMoC>6 shows multiple features, indicating an 



unusual electronic state of Mo in this compound. In or- 
der to explore whether the multiple feature within the 
main peak region of Mo 3d spectrum is a consequence of 
a negative U e f / , we have performed a core level calcula- 
tion considering a cluster with four Mo atoms arranged 
at the corners of a square. Each atom had 4 orbitals, with 
a crystal field splitting of 0.9 eV between a triply degen- 
erate set (t2 S -hke) and a single e 9 -like orbital. Hopping 
between ti g orbitals was taken to be -0.4 eV, while that 
between the e g -like orbitals was taken to be -0.5 eV. U e ff 
was taken to be -0.5 eV. The presence of the core- hole was 
taken into account by a decrease in the onsite energies 
by the core-hole potential of 0.9 eV The resulting 

calculated spectrum is shown in the main figure as a solid 
line for comparison. The agreement between the experi- 
mental and the calculated spectra is remarkably good. It 
was also found that the calculated spectrum had invari- 
ably the simple double peak structure of Mo and M0O3 
for any parameter set as long as the U e f / was taken to be 
positive, indicating that the U e ff at the Mo sites in in- 
deed negative in this system, driving an electronic phase 
separation. 

Combining ab-initio and model Hamiltonian calcula- 
tions, we establish the presence of strongly enhanced J e / / 
and J, responsible for the unusual electronic and mag- 
netic properties of S^FeMoOe- The enhanced / results 
in an effective negative U at the Mo site which results in 
a complex core level spectrum observed in our photoe- 
mission experiments. 
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I. FIGURE CAPTIONS 

Fig. 1. LMTO band dispersions for Si^FeMoOg along 
symmetry directions (see text). 

Fig. 2. The dispersions of the bands with t2 g symmetry 
at the Mo site for (a) up spin and (b) down spin. All other 
degrees of freedom have been integrated out as discussed 
in the text. 

Fig. 3. Energy level diagram expected at the Fe (left 
panel) and Mo (central panel) sites from an ionic picture. 
The energy levels at Mo are modified in the presence of 
Fe-Mo hopping interactions as shown in the right panel. 

Fig. 4 The experimental (solid circles) and calculated 
(solid line) Mo 3d core level spectrum for S^FeMoOV 
The experimental spectra (solid circles) for for M0O3 and 
Mo-metal are also shown. 
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